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A unidimensional model is developed including the equations of continuity and force equilibrium and a physi¬ 
cal equation describing the state of glass melt in the fiber formation zone. Perturbations related to variations in 
the technological parameters are introduced into the solutions of equations detennining a stable process, and 
possible deviations of the fiber radius are found. 


Stringent requirements are imposed on optical fiber, in¬ 
cluding stability of fiber diameter. 

Analytical models constructed in the context of fluid dy¬ 
namics make it possible to predict the effect of variation of 
technological parameters (the rate and force of drawing, the 
temperature and viscosity in the deformation zone etc.) on 
stability of the drawing process. 

The scheme of the process considered is shown in Fig. 1. 
Fiber 1 is drawn from glass melt 2 heated using a particular 
method (in a furnace, with a laser, etc.). According to the 
heating method, glass can exist as an intermediate product or 
as a melt flowing out of the opening in the crucible bottom. 

The glass deformation zone can be divided into two 
zones: molding zone a and drawing zone b. 

The molding zone is characterized by a fast narrowing 
and then it gradually passes into the drawing zone. In this 
zone the fiber characteristics are finally formed, which are 
influenced by perturbations both from zone a and from the 
drawing mechanisms beneath, for instance, in the form of 
mechanical vibrations. 

Zone b can be studied using unidimensional models, in¬ 
cluding continuity and force equilibrium equations and a 
physical equation describing the state of the glass melt [1], 
Without taking into account elastic properties of the glass 
melt, these equations have the following form [2]: 
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where r(z, t) is the glass melt flow radius; t is the time; z is 
the axial coordinate; v (z, t ) is the glass melt flow velocity; 
p is the glass melt density taken to be constant; p (z, t) is the 
stress in the considered flow section; p is the dynamic vis¬ 
cosity. 

Substituting expression (3) in condition (2), we obtain 

djvr 1 ) | d(v 2 r 2 ) ^ d(\xr 2 dv / dz) _ Q 

dt dz dz 

The stationary variants (1) and (4) have the following 
form: 

r 2 v=Q = const; (5) 

3pr 2 — - pvV 2 = C, (6) 

dz 

where C is the integration constant. 

Without taking into account the inertia term, Eq. (6) can 
be represented using expression (5) in the form 

dv 

— = Ndz, (7) 
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Fig. 1 . Glass melt deformation zone. 
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where N= --— ; p 0 and v 0 are the radius and the flow ve- 

3(ir 0 2 v 0 

locity, respectively, at z = 0. 

Equations (5) and (7) make it possible to determine the 
variations of the flow velocity along the drawing zone and 
the shape of this zone. The length of the drawing zone is L. 
At z = L the flow velocity is equal to the drawing velocity v d , 
which is generated by the drawing mechanism. 

Let us assume that in the drawing zone p = const. Then 
from Eqs. (5) and (7) we get 

v = v 0 explain k); (8) 

r = r 0 ex (9) 

where £, = z/L; k = v d /v 0 is the neck coefficient. 

The solutions obtained describe a stable drawing process. 
To estimate the effect of variations in the technological 
parameters on the stability of the drawing process, let us in¬ 
troduce perturbations of the first order for the radius r, ve¬ 
locity v, and viscosity p, and let us seek the first-order solu¬ 
tions of Eqs. (1) and (4) in the form 

r (\) =r(z)\\ + 7(z,0]; 

< v (1) = v(z)[l + v(z,t)]; (10) 

b(i) =P(")[ 1 + P(zd)], 


where r (z ), v (z ), and p (z) are solutions determining a sta¬ 
ble process. 

Substituting expression (10) in condition (1) and neglect¬ 
ing terms of the second order of smallness, we obtain 
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Without considering the inertia term, Eq. (4) after substi¬ 
tuting conditions (10) into it will take the form 

v z —+ v+2 r + v = P(t), (12) 

v z 


where v, = — ; v z = — ; v = — is the kinematic viscosity; 

dz dz p 

P (t ) is a certain arbitrary time function. 

Equations (11) and (12) can be regarded as perturbation 
equations. The solution of these equations depends on the 
problem setting. For example, let us set the viscosity pertur¬ 
bation v caused by temperature fluctuations and let us find 
the respective perturbations r of the fiber radius. For this 
purpose we substitute solutions (8) and (9) for a stable pro¬ 
cess into Eqs. (11) and (12). As a consequence we obtain 

7" + 7'e“ 2a = v', (13) 


where r" 
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Fig. 2. The effect of viscosity on the deviation of normalized fiber 
radius. 


Solutions of Eq. (13) depend on the boundary conditions 
and on the form of the function determining the viscosity 
perturbation. 

Let us assume that the conditions of feeding glass melt 
into the drawing zone and gripping the drawn fiber do not 
change. Then one can take 

7 (0) = 7 (a L ) = 0, (14) 

where a L is a parameter correlating with q = 1. 

Let us choose the exponent as the function determining 
the viscosity variation. Let us assume that there is no pertur¬ 
bation at a = 0. Viscosity perturbation is determined by the 
function 
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where v 0 and P are parameters determining the disturbance 
amplitude. 

Let us consider the stationary variant of Eq. (13). In this 
case the perturbation r nonnalized with respect to v 0 for 
function (15) will take the form 

^-IgWl + J+Cjtx + Cj, ( 16 ) 
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where C and C 2 are integration constants determined by 
boundary conditions (14): 


C2= F 
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where b = In (k/2). 

Figure 2 shows the results of calculation of the norma¬ 
lized disturbance of radius F/v 0 based on expressions (16) 
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and (17) for different values of p. Let us correlate the maxi¬ 
mum disturbance of viscosity v max and radius r max and de¬ 
rive a relationship based on the calculation results: 
Vmax/Lnax=” : P = 1 , n = 6; P = 2 > « = 9.2; P = 4, n = 4.6; 
P = 6, n = 3.8. 

In the most favorable case (P = 2) a 10% deviation of 
viscosity from a preset value will produce an approximately 
1% deviation from the preset fiber radius, and in the worst 
case (P = 6) an approximately 3% deviation of the radius. 


Thus, using Eqs. (11) - (13), it is possible to estimate the 
effect of the perturbation of the drawing velocity v on the de¬ 
viation of the fiber radius r. 
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